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1. 


INTRODUCTION 


Hicu polymers generally exhibit a discontinuity in thermal expansion at a 
characteristic temperature. This temperature also marks a relatively abrupt 
change in other thermal, mechanical, electrical and optical properties. 
Ehrenfest! has termed these discontinuities as second order transitions. A 
considerable body of experimental data on high polymers has been reported 
by a number of investigators in this direction. Recently Subrahmanyam,” 
Krishnamurthi and Sivarama Sastry,’ during their study of the variation of 
ultrasonic velocities with temperature in polymerized methyl methacrylate 
and polystyrene, have reported from this laboratory similar transitions. 
In the present investigation, these two polymers have been studied by photo- 
elastic methods. 


2. DESCRIPTION OF THE APPARATUS AND METHOD OF OBSERVATION 


All the specimens needed were cut from cast sheets of polymethyl metha- 
crylate or perspex and polystyrene. The cut surfaces of the samples were 
filed flat and smooth, and the samples were annealed at 90° C. for about 40 
hours. The specimens were then worked out as plane parallel plates and 
were well polished. None of the specimens showed double refraction under 
crossed nicols in an unstrained condition. 


A parallel beam of light from a sodium-vapour lamp was passed through 
a Nicol, the vibration direction of which was maintained at 45° to the vertical. 
The emergent beam then passed through the test-piece into a Babinet compen- 
sator, the principal axes of which were vertical and horizontal. After 
passing through the compensator the light was observed through an eyepiece 
containing a Nicol crossed with respect to the polarizer, and the usual 
Babinet fringes were obtained. When the specimen was compressed, the 
fringes shifted one way or the other and the magnitude of the shift was a 
Measure of the path difference produced between the horizontally and verti- 
cally vibrating beams of light passing through the compressed specimen. 
Compression was produced by a lever arrangement. 
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Measurements of the fringe shift were made for stresses of about 
100 kg./cm.? in the case of perspex and 60 kg./cm.” in the case of poly- 
styrene. Following Pockels, the displacement of the fringes was measured at the 
middle of the length of the prisms at three places along its breadth, namely, 
near the left edge, the middle, and the right edge. The final measurements 
were made after ascertaining that the distribution of load was fairly uniform. 
The mean shift of the Babinet fringe in each case was obtained by taking the 
average of several readings. The path retardation read in terms of the head 
scale divisions of the Babinet compensator is a measure of the birefringence 
produced in the samples. 


In order to study the variation of the induced birefringence with tempe- 
rature, a compact electrical heater was built, the specimen being centrally 
situated within the enclosure. By means of a variac the voltage could be 
varied slowly and steadily till the requisite temperature was attained and a 
measurement on the Babinet was made only after ascertaining that the tempe- 
rature was constant for about two hours. The procedure was repeated till 
consistent results were obtained at each temperature both while heating and ~ 
cooling the specimens. Measurements of fringe shifts were made at intervals 
of 10° C., temperatures being read by a sensitive thermometer. The mean 
values of the Babinet shifts so obtained at each temperature were then made 
use of in deducing the relative stress-optical coefficients. 


As the refractive index enters in the evaluation of the photoelastic 
constants, it is found necessary to determine the same for perspex and 
polystyrene at various temperatures. In the present investigation the usual 
method‘ in which dn/dt is directly measured as a function of temperature 
is employed. The method consists in cutting a plane parallel plate of thick- 
ness / out of the specimen and obtaining localized interference fringes by 
reflection from the two surfaces of the plate at normal incidence. The 
specimen was kept in a suitable furnace and the fringes, which consisted of 
sharp, bright and almost circular rings, were observed by means of a micro- 
scope. When heated, the fringes began to expand in both the cases and a 
continuous observation of noting the temperature corresponding to the 
passage of successive fringes was carried out. The rate of heating was so 
slow that it took eight hours to heat the specimen from room temperature to 
90°C. The observations were repeated with cooling specimen, this also 
taking the same time. 


The measurement of the expansion gives one the alteration in the path. 
Knowing At corresponding to each wavelength shift, dn/dt was calculated 
for the mean temperature making use of the formula 
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dn/dt = 41. dx/dt — na, 


where a is the coefficient of thermal expansion at temperature t. The values 
of a at various temperatures were available from the work of Yasaku Wada. 
The absolute values of indices of refraction (mp) for sodium light at room 
temperature have been determined to three significant figures using the material 
in the form of tiny prisms, employing a spectrometer. In both the cases 
prisms having an angle roughly of about 60° were cut from sheets, 6-5 mm. 
in thickness, and the two sides of the prisms making the angle were well 
polished. Making use of these values of (np) at 25° C., the refractive indices 
at various temperatures were calculated from the data for dn/dt by numerical 
integration. 
TABLE I 
Values of Xp at different temperatures 


Temperature ° C. Perspex Polystyrene 
25 1-49300 1-59500 
30 1-49047 1-59328 
40 1-48542 1-58986 
50 1-48043 1-58646 
60 1-47539 1 -58308 
70 1-45906 1 -57973 
80 1-45193 1 -54927 
90 1-44447 1-54113 


The calculated values of the refractive indices (mp) at various temperatures 
for perspex and polystyrene are given in Table I. 


3. RESULTS AND DISCUSSION 


The photo-elastic behaviour of isotropic substances is given by the 
equations : 


By — B= —q-ZZ; -B=—4_-ZZ; — B= — 4y,-ZZ; and 
By = Bs, = B,, = 0. 

In these equations B = 1/n®, B,,? = 1/n?,,, etc., where n is the refractive 
index of the undeformed specimen of an isotropic substance, 1), etc., are the 
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components describing the Fresnel ellipsoid in the deformed condition, and 
ZZ, a simple linear compressional stress along the z-axis. 4; and qj. are 
the only two independent stress-optical coefficients, (411 — 412) being equal 


The determination of the relative stress-optical coefficients (q,, — qs) 
was carried out by employing a Babinet compensator, the observations being 
given in detail in Tables II and III. 470-5 divisions of the Babinet compen- 
sator corresponded to a path retardation of one wave-length of sodium light. 
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TABLE II 


Results for perspex 


Mean shift in (Gir — 412) (Pir — Pre) 
Temperature °C. divisions of x 1078 x 104 
head-scale 

27 104-0 2-626 1-133 

40 110-0 2-846 

51 116-0 3-031 

61 121-0 3-271 

71 130-0 3-572 

80 139-0 3-922 

91 150-0 4-321 

TABLE III 
Results for polystyrene 
Mean shift in — 412) (Cr — C12) (Pir — Pra) 
Temp. ° C. divisions of x 104% xi x 104 
head-scale 

32 144-0 4-733 2-782 1-371 
42 137-0 4-540 2-722 1-236 
52 130-0 4-347 2-672 1-161 
62 122-0 4-115 2:614 1-075 
72 113-0 3-903 2-552 0-996 
81 100-0 3-678 
92 84-0 3-157 
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From the mean shift of the Babinet fringe, the path differences corres- 
ponding to unit stress (1 dyne cm.~*) and unit length (1 cm.) of light beam in 
the stressed specimen were calculated knowing that the path difference is 
5-893-10->cm. for a shift of 470-5 divisions. The stress-optical constant 
in the corresponding direction of observation is 25/n®. In the case of poly- 
styrene the shift was in the same direction as for common glass while it was 
opposite in the case of perspex. Hence, the sign of (41: — qi) is negative 
for polystyrene and positive for perspex. Also an examination of Tables 
II and III reveals that the numerical value of the stress-optical constant 
(41. — 412) for perspex increases with increase of temperature while it 
decreases in the case of polystyrene. Furthermore the values of the stress- 
optical constant at the lowest temperature differ by 60 to 70 per cent. from 
those at about 90°C. 


The values of the strain-optical constants (p,; — P12) are calculated from 
the equation connecting p’s, q’s and the elastic constants, namely, 


(Pir — = (C1 — Cra) (Gr — 912). 


The elasticity data on polystyrene at different temperatures have been 
obtained by Subrahmanyam.? The rigidity modulus of perspex is not 
available in the literature. Robinson, Ruggy and Elizabeth Slantz* have 
studied the stress-optical properties of polymethyl methacrylate as a func- 
tion of temperature. However, their results do not give the complete data 
of the photo-elastic constants. 


Temperature variation of the birefringence produced by stress in the 
above plastics is represented graphically in Fig. 1. A pronounced transi- 
tion in the birefringence vs. temperature graphs of perspex and polystyrene 
will be noticed at 62°C. and 70°C. respectively. As has already been 
mentioned Subrahmanyam,” Krishnamurthi and Sastry* of this laboratory 
have studied the variation of ultrasonic velocities in perspex and polystyrene 
with temperature. They observed that both perspex and polystyrene show 
a transition at 61°C. and 70°C. respectively. The samples used by the 
authors and the above workers are taken from the same sheet of material. 
From the above facts it can be concluded that the present photo-elastic 
investigations in these two thermoplastics clearly show the existence of the 
transition points in good agreement with the earlier workers. 
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4. SUMMARY 


The birefringence produced in samples of polymethyl methacrylate 
and polystyrene when subjected to stress is measured at different tempera- 
tures employing a Babinet compensator. Graphs of birefringence against 
temperature are plotted for the temperature range 30 to 90° C. for perspex 
under 100 kg./cm.? stress and for polystyrene under 60 kg./cm.? The bi- 
refringence of the stressed methacrylate is positive and for polystyrene it is 
negative. The birefringence in the case of perspex increases with temperature 
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while in polystyrene it shows a decrease as the temperature increases. There 
is a transition point at 62° C. for perspex and at 70° C. for polystyrene. The 
relative stress-optic and strain-optic constants are evaluated at each tempe- 
rature taking into consideration the temperature variation of the refractive 
index. 
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(Communicated by Sir C. V. Raman) 


1. INTRODUCTION 


ACCORDING to Dana,! magnesite (MgCO,) commonly occurs in massive and 
compact forms but clear single crystals are rare. It is not surprising there- 
fore that few investigations have been made on its Raman spectrum. The 
only study reported in the literature was by Schefer and others* who, using 
a single crystal, have reported a line at a frequency shift of 1096cm.-! The 
opportunity for the present study arose due to the fact that some clear crystals 
from Brumadas, Bahia, Brazil were made available to us through the kind- 
ness of Dr. C. S. Pichamuthu. The density of a clear specimen was deter- 


mined and found to be 3-02 gm./c.c., in fair agreement with the value 3-037 
given in the International Critical Tables. 


As the specimens were colourless, it was thought that the complete 
Raman spectrum could be readily recorded using the A 2536-5 resonance line 
of mercury for excitation. Actually however, the crystals were found to be 
opaque to radiations of shorter wave-length than A 2652. Hence the visible 
radiations of the mercury arc and a glass spectrograph had to be used, and 
the recording of the complete Raman spectrum was found to be far from 
easy, since in the first place the crystals were of rather small size and not 
completely free from flaws and since they also exhibited a faint luminescence. 
A strong and sharp line at a shift of 332 cm. and a strong but slightly diffuse 
line at 1096 cm.-! with a width of about ten wave-numbers were the principal 
features recorded on all the plates. By interposing a filter of NaNO, solu- 
tion to suppress A 4046 from the incident radiation and by placing a thin 
opaque strip just before the photographic plate to screen the exciting radia- 
tion A 4358, clearer spectra were recorded in which could be observed and 
measured a line with a frequency shift of 212cm.—! In some particularly 
clear spectrograms, a weak and rather diffuse line was found with a frequency 
shift of 735cm.-! By analogy with the Raman effect in calcite one should 
expect one more line with a frequency shift in the vicinity of 1450cm-! A 
faint trace of such a line was found on one of the plates with an estimated 
frequency shift of 1460cm.-? 
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2. DISCUSSION OF RESULTS 


Magnesite has the same crystal structure as calcite belonging to the 
space-group D%, of the trigonal system. Dolomite has also a similar struc- 
ture, though it belongs to the space-group C2;, In view of the similarity 
of crystal structure and chemical composition it is of interest to compare 
their Raman spectra. Table I shows the observed frequency shifts and their 
relative intensities in the case of calcite, dolomite and magnesite. The data 
for calcite are from a paper by Krishnan® while those for dolomite are due 
to Couture by whom the Raman effect in dolomite has been investigated in 
detail. 
TABLE I 
Raman Frequencies 


Substance Frequency shifts in cm.-! 


284 (15) 


CaCO, 156 (8) 712(2) 1086 (20) 1434(5) 
CaMg (COs)s| 176(9) 301 (10) 335 (0-3) 724(2) 1099(23) 1444 (1-4) 
MgCO, ..| 212(w) 332(st) 735(w) 1096 (st)  1460(2) 


As is to be expected, Table I exhibits a close similarity between the Raman 
spectra of calcite, dolomite and magnesite. All the frequencies of magnesite 
are definitely higher than those observed with calcite, this feature being more 
pronounced for the lattice oscillations than for the internal frequencies of 
the CO; group. An explanation of this feature is to be sought for in the fact 
that the lattice constants of calcite and magnesite are respectively 6-36 A.U., 
and 5-61 A.U., indicating a notable decrease in the interionic distances in 
the latter case. This would result in an increase of the interionic forces and 
hence also of the vibration frequencies involving the CO, groups. 


Schefer et al.® have studied the infra-red reflection and absorption 
spectra of magnesite and observed the following frequencies amongst others: 
Reflection Spectrum: 726, 889 and 1495cm.-!; Absorption Spectrum: 
730, 877, 899 and 1430cm.-! The only study on the far infra-red reflection 
spectrum of magnesite was by Coblentz® who observed a reflection maximum 
at 30-7 p, i.e., 326cm.- It will be noticed from the data that the frequency 
observed at 889 cm.-! in the infra-red does not appear in the Raman effect, 
whereas the frequency 1096cm.-! observed in the Raman effect does not 
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appear in infra-red absorption. These facts are in accordance with the 


theoretical analysis of the vibrations of the calcite structure by Bhagavantam 
and Venkatarayudu.’ 


The author wishes to express his sincere thanks to Prof. Sir C. V. Raman, 
F.R.S., N.L., for suggesting this investigation and for the kind interest he 
took in it. 

SUMMARY 


The Raman spectrum of magnesite (MgCO,) has been investigated and 
found to be analogous in its features to that of calcite, as is to be expected 
in view of the similarity of structure. The frequency shifts observed are: 
212 (w), 332 (st), 735(w), 1096 (st) and 1460(?)cm.-! All the frequencies 
and especially the two lowest are notably higher than for calcite and this is 
explicable as due to the smaller interionic distances and correspondingly 
larger force constants. 
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Tue reason for the marked instability of ceric halides have not been fully 
understood. Only the fluorides are known as normal salts of both cerous 
and ceric cerium.!* Ceric chloride has not been isolated in the free state. 
When ceric hydroxide is dissolved in concentrated hydrochloric acid, a dark- 
red solution containing ceric chloride or the acid H,CeCl, is obtained which 
rapidly decomposes with evolution of chlorine. The possibility of obtaining 
ceric chloride or hexachloroceric acid in absence of water was first reported 
by Koppel* and following his method, several double salts of the acid with 
organic bases have been prepared.*’ It is of interest to note that while the 
sulphato- and nitrato-ceric acids are stable in aqueous solutions, the chlorato- 
ceric acid is quite unstable. Although Koppel, for the first time, tried to 
isolate ceric chloride by carrying out the reaction between ceric oxide and 
hydrogen chloride in methyl alcohol, he only succeeded in preparing the 
pyridinium double salt of hexachloroceric acid. This paper deals with a 
systematic study for the isolation of hexachloroceric acid by using methyl 
alcohol and dioxane as reaction media to increase the stability of the products. 


The reaction of hydrated ceric oxide and hydrogen chloride in methyl 
alcohol did not yield useful results in the isolation of hexachloroceric acid. 
The failure to isolate the reaction product from methyl alcohol solution 
could be either due to the high solubility of the product in methyl alcohol or 
due to the incapacity of the solvent to co-ordinate with the product as a 
solvate to yield a stable complex or due to both. Dioxane was next employed 
as the solvent. Dioxane was expected to be a good solvent to prepare ceric 
chloride since the solvent might increase the stability of the compound by 
forming a type of oxonium compound. Several inorganic halide dioxinates 
have been reported.*!* 


It has been found in the present investigation that hexachloroceric acid 
crystallises out, under suitable conditions from a dioxane solution with four 
molecules of dioxane for every mole of the acid. The exact nature of the 
mechanism of reaction of hydrated ceric oxide and hydrogen chloride in dio- 
Xane and the formation of the crystalline product is not understood with 
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certainty. Although the formation of hexachloroceric acid is noticed when 
various other organic solvents are used as reaction media, viz., methyl alcohol, 
ether, ethylene glycol mono methyl ether, etc., the crystalline product is 
obtained only in presence of dioxane. Evidence has been adduced to suggest 
that an anion complex of the type CeCl, is absent but that the product is a 
molecular compound of ceric chloride and hydrogen chloride of the formula 
CeCl,. 2 HCl, co-ordinated with four molecules of dioxane. If the hexa- 
chloroceric acid and its pyridinium salt were complexes with cerium attached 
to the six chlorine atoms in a stable anion, silver nitrate solution (in non- 
aqueous medium) should give no precipitate. Again, if it is assumed that the 
anion consists of all the chlorine atoms combined with the cerium atom, 
there is no reason why one chlorine atom should behave differently in liberat- 
ing iodine from potassium iodide. Studies on the physico-chemical proper- 
ties like conductivity, molecular weight, vapour pressure and _ spectral 
characteristics—details of which are to be published in another part of the 
series—strongly support the view that the compound does not contain an 
anion complex of the type CeCl,.-- 


EXPERIMENTAL 
Materials used 


(1) Hydrated ceric oxide.—Pure ceria was prepared from ceric ammonium 
nitrate by the iodate method. The ceric oxide thus obtained was converted 
into nitrate, dissolved in water and made alkaline with dilute ammonia. The 
precipitated hydroxide was oxidised with hydrogen peroxide and let stand 
for 12 hours. The yellow product was filtered, washed thoroughly with 
water, dioxane and petroleum ether when a free flowing yellow powder 
(CeO, -2 H,O) was obtained. 


(2) Hydrogen chloride——Hydrogen chloride gas was prepared by the 
method of Harwood Tucker.** 

(3) Methyl alcohol, dioxane and petroleum ether (30-75°C.) were 
purified and dried by standard methods.™ 
Analytical 


(1) Cerium.—Cerium was precipitated as the oxalate, ignited and weighed 
as ceria. 


(2) Active chlorine—A known quantity of the material was put into 
acidified potassium iodide solution and the iodine liberated was estimated 
by standard sodium thiosulphate. 


(3) Total chlorine—A known weight of the substance was dissolved in 
1 N sodium hydroxide solution and the contents acidified with chloride-free 
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nitric acid until the precipitated cerous hydroxide redissolved. The decom- 
position of the hypochlorite was effected by warming the solution with a 
little urea. The chloride content of the solution was estimated volumetri- 
cally by the silver nitrate method. | 


(1) Reaction of hydrated ceric oxide and hydrogen chloride in methyl alcohol 


Hydrated ceric oxide (5g.) and methyl alcohol (100c.c.) were taken 
in a flask at 20° C. and hydrogen chloride was passed for nearly three hours 
with constant stirring, at the rate of about 6 litres per hour. The methyl 
alcohol solution turned orange-red in colour. The flask was kept at 2-3° C. 
for two days when the solution became more intensely coloured. The solu- 
tion was filtered, out of contact with moisture and the residue on analysis 
was found to consist of 97-:1% CeO, .2 H,O and 3% CeCl,. The filtrate 
was analysed for active chlorine and cerium contents. The estimation of 
total chloride was not attempted as it was not possible to free the solution 
completely of dissolved hydrogen chloride. It was noticed that for every 
gram atom of active chlorine, a corresponding quantity of cerium was present 
in the methyl alcohol solution. These results in general indicated that about 
3% of ceric oxide was converted into ceric chloride or hexachloroceric acid 
which dissolved in alcohol producing the orange-red colour. 


The above experiment was repeated at 0°C., — 15°C. and — 75°C. 
In the first two cases, there was not much difference in the rate of reaction 
and development of colour from the previous experiment. There was 
practically no reaction at — 75°C. The reaction could not be carried out 
for a long time as the methyl alcohol solution became very viscous, prevent- 
ing the free passage of hydrogen chloride. 


The product obtained by the reaction of 7-5g. of ceric oxide and 
hydrogen chloride in 150c.c. of methyl alcohol for three hours, was kept 
at 0°C. The progress of reaction was followed by estimating the active 
chlorine in aliquots of the solution every twelve hours. From these values, 
the amount of cerium (IV) in the solution was calculated. The results are 
presented in Table I. 


These results show that in thirty-six hours, the ceruim (IV) content has 
increased nearly seven-fold, after which the rate of reaction is very slow. 
Attempts at the isolation of hexachloroceric acid from methyl alcohol solution 


About 100 c.c. of the ceric chloride in methyl alcohol, was subjected to 
vacuum desiccation when the intensity of colour and viscosity of the alcohol 
solution were found to increase. During the final stages an yellowish-red 
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solid separated which on further desiccation yielded only cerous chloride. 
Presence of free chlorine in the exit vapour was detected towards the final 
stages of evaporation. Concentration by cooling at 0° C. and at — 20°C. 
gave no better results. 


TABLE I 


Progress of reaction between ceric oxide and hydrogen chloride in 
methyl alcohol with time 


Time % Ce (IV) 
(in hours) in soln. 


0-33 
1-25 
1-89 
2°22 
2°36 
2°58 


72 2-68 


Removal of methyl alcohol from the solution was also attempted by 
passing a current of hydrogen chloride to maintain the acidity at a high value. 
The solid that separated was found to be merely cerous chloride. Separa- 
tion of the ceric chloride by fractional freezing of the alcohol did not result 
in higher concentration of the ceric chloride. Fractional precipitation of 
the chloride by the addition of benzene, carbon tetrachloride, chloroform, 
ether, dioxane and petroleum ether did not yield any fruitful results. It can 
therefore be concluded that methyl alcohol is not a suitable medium for the 
preparation of hexachloroceric acid. 


(2) Reaction in presence of dioxane 


Hydrated ceric oxide (5 g.) suspended in dioxane (100c.c.) was next 
treated with hydrogen chloride for nearly three hours. As there was consi- 
derable rise in temperature, the flask was cooled in water at 25°C. The 
development of the orange colour in the solution was much quicker than 
with methyl alcohol and left only a small quantity of residual solid. The 
flask was stoppered and left overnight at 2-3°C. when some orange-red, 
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needle-shaped crystals were formed. At the same time the colour of the 
solution had become paler. The solution was filtered through a sintered 
glass funnel, out of contact with moisture. The filtrate on analysis was 
found to contain cerium and active chlorine in the gram atom ratio of 1: 1. 


The residue in the funnel consisting of the crystalline mass and unreacted 
ceric oxide was dried in a current of dry air after washing with a small quantity 
of dioxane to remove the adsorbed solution and hydrogen chloride. The 
solid was then extracted with 50c.c. of dioxane and filtered when a deep 
red filtrate was obtained, leaving a white residue on the filter. The white 
residue thus obtained was found on analysis to be a mixture of 87% ceric 
oxide and 13% CeCl,. The deep red filtrate was then analysed for cerium, 
active chlorine and total chlorine and the ratio was found to be 
1:0:994:6-10. From this, it was concluded that the crystals had the 
empirical formula H,CeCl, with or without dioxane of crystallisation. 


In subsequent experiments, 10 g. of ceric oxide were treated with hydro- 
gen chloride, using the dioxane solution left as the filtrate, andthe unreacted 
ceric oxide was removed by filtration. The filtrate separated into two layers, 
the heavier one being deep red and the lighter one orange-red in colour. 
On cooling the heavier layer to 2-3° C. for about 6 hours, a large crop of 
needle-shaped, bright, orange-coloured crystals separated, leaving behind 
pale orange coloured mother liquor. It was also noticed that the presence 
of hydrogen chloride was necessary for efficient crystallisation. 


Analysis of the crystals 


The crystals prepared as above were washed several times with petroleum 
ether and dried in a current of dry air, and analysed for cerium, active chlo- 
tine and total chlorine. The mean values of six different crops of crystals 
were: Ce = 19-48%; active chlorine = 4:94%; and total chlorine = 
29 69%. Calculated percentage composition of H,CeCl,-4 C,H,O, is: 
Ce = 19-80%; active chlorine = 5-02% and total chlorine = 30-13%. 
The cerous chloride obtained by thermal decomposition of the crystals was 
also analysed. A known weight of the crystals was heated to 150°C., out 
of contact with moisture and cooled in a current of nitrogen and weighed. 
It was then analysed for cerium and chloride. The average value was 
33-62% residue analysing to 56-18% Ce and 42-85% Cl. Calculated value 
of CeCl, from HyCeCl., 4 CyH,O, is 34-87 % with 56-79% Ce and 43-21% Cl. 
These results conclusively show that the residue on decomposition of the 
acid crystals is cerous chloride, and the per cent. weight of CeCl, obtained 
agrees with the calculated value when the acid has the formula 
H,CeCl,-4 C,H;O,. Because of the high reactivity of the crystals towards 
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light and moisture, experiments to estimate carbon and hydrogen did not 
give consistent results. 


Preparation of salts 


Koppel prepared the organic derivatives of hexachloroceric acid by 
treating the methyl alcohol solution with organic bases. The hexachloro- 
ceric acid crystals as prepared above, were also employed to prepare these 
compounds as follows. To a solution of the crystals in dioxane was added 
a solution of pyridinium chloride in methyl alcohol. Yellow shining crystals 
were thrown out, the amount of which increased by the addition of ether. 
The crystals were filtered off, washed with ether and analysed. The per- 
centage composition was Ce = 26-80, active chlorine = 6-72, total 
chlorine = 39-26 and nitrogen = 5-12 agreeing well with the formula 
(C;H;NH), CeCl,, for the complex. 


With a view to prepare the sodium salt, the following methods were 
tried. Samples of methyl alcohol solution of the acid crystals were 
added to a methyl alcohol solution of sodium hydroxide, a solution of sodium 
methoxide in methyl alcohol and dry solid sodium chloride respectively. 
In the first two cases, there was immediate decolourisation with precipitation 
of cerous chloride. In the experiment where sodium chloride was added 
there was no immediate change, but on standing the solution got slowly 
decolourised, with deposition of white residue. Further studies on the 
preparation of salts of the acid are in progress. 


Attempts to prepare the acid without dioxane of solvation 


A portion of the heavy layer of liquid that is formed during the reaction 
of ceric oxide and hydrogen chloride in dioxane at 50° C. was subjected to 
evacuation at 25°C. and 0-1 mm. pressure. As the solid dried up, 
a pale yellow powder was left which was identified as cerous chloride. 
Repeated washing of the liquid with petroleum ether also resulted in 
getting cerous chloride as the residue. When the acid crystals were 
subjected to evacuation, decomposition took place leaving behind cerous 
chloride. Gentle heating of the crystals also resulted in decomposition. It 
was thus concluded that the acid got easily decomposed when the 
dioxane of solvation was removed. 


SUMMARY 


(1) The reaction of hydrated ceric oxide with hydrogen chloride in 
methyl alcohol under controlled conditions did not result in the isolation of 
hexachloroceric acid. 


Studies on Hexachloroceric Acid—I 219 


(2) Hexachloroceric acid has been obtained as a crystalline product 
from the reaction of hydrated ceric oxide and hydrogen chloride in dioxane. 


(3) Conditions have been worked out for obtaining the crystals in a 
pure form. Analytical results agree with the formula H,CeCle, 4 C4HsOx. 
It has not been possible to obtain the acid without the dioxane of solvation. 


(4) The pyridinium salt of the acid has been prepared, the analytical 
values agree with the formula (C;H;NH),CeCl, for the compiex. 


The authors’ thanks are due to Prof. K. R. Krishnaswami for his keen 
interest and helpful suggestions during the course of the work. One of the 
authors (S. S. Moosath) is indebted to the University of Travancore for 
deputation to the Indian Institute of Science, Bangalore, to conduct this 
work. 
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STUDIES ON HEXACHLOROCERIC ACID 


Ii. Reaction of Hydrated Ceric Oxide with Hydrogen Chloride Gas in 
Non-Aqueous Solvents 


By S. S. MoosaTH 
(Department of General Chemistry, Indian Institute of Science, Bangalore) 
Received February 25, 1956 


(Communicated by Dr. M. R. A. Rao, F.A.Sc.) 


RESULTs of the reaction of hydrated ceric oxide and hydrogen chloride in 
methyl alcohol and dioxane have been detailed in Part I of the series. The 
reaction was studied in various liquids—both polar and non-polar—which 
are usually good solvents. Reaction was found to take place in some of 
the solvents, although in the absence of dioxane, no crystalline product 
could be isolated. The reaction is found to take place in such solvents 
which can form oxonium type of compounds with hydrogen chloride. The 
failure to isolate hexachloroceric acid from solutions in methyl alcohol, 
ether and ethylene glycol mono-methyl ether, in the absence of dioxane 
clearly shows that a suitable solvate is necessary for the complex formation 
in a crystalline form. Addition of dioxane to concentrated solutions in 
various solvents results in the formation of crystals having the composition 
H,CeCl,, 4C,H,O.. It appears that the solubility of the product is very 
much enhanced by the presence of OH group in the solvent as evident from 
the high solubility in methyl alcohol. The extremely rapid reaction in 
ethylene glycol mono-methyl ether and the non-separation of a second layer 
or crystals, support the view that the presence of tetravalent oxygen favours 
the reaction and that the presence of OH group in the solvent increases the 
solubility of the product formed. The separation of a second layer in ether 
shows the relatively low solubility of the product in ether which does not 
have the hydroxyl group. The exact nature of the mechanism of reaction 
in these solvents is not understood with certainty. 


EXPERIMENTAL 


Material used.—Hydrated ceric oxide and hydrogen chloride were pre- 
pared by the method detailed in Part I of the series. Dioxane, petroleum 
ether, benzene, chloroform, acetone, ether, ethylene glycol mono-methyl 
ether, formamide, and cyclohexanol were purified and dried by the standard 
methods.?» 2 
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Reaction in acetone, cyclohexonal, formamide, anisole benzene, and carbon 
tetrachloride : 


(a) Acetone.—A mixture of hydrated ceric oxide (5g.) and acetone 
(100 c.c.) was subjected to the action of hydrogen chloride for three hours. 
The acetone turned only light yellow in colour. When the flask was left at 
2-3° C. for a day, the liquid became dark brown in colour with the forma- 
tion of mesitylene oxide due to the polymerisation of acetone in presence of 
hydrogen chloride. Acetone was therefore considered unsuitable for the 
preparation of hexachloroceric acid. 


(b) Cyclohexanol.—The previous experiment was repeated with cyclo- 
hexanol as the reaction medium. The ceric oxide remained unreacted. 
The solvent became slightly brownish in colour but no crystals separated on 
cooling the solution. 


(c) Formamide.—When this solvent was used as the reaction medium, 
there was no change in the ceric oxide. The formamide became slightly 
brownish in colour but very little of the hydrochloroceric acid was produced. 


(d) Anisole-—In this case also, the ceric oxide remained practically 
unreacted. But on keeping, the liquid became dark brown but no hydro- 
chloroceric acid was produced. 


(e) Reaction in benzene and carbon tetrachloride—These solvents were 
tried as reaction media. As before, hydrated ceric oxide was suspended 
in these solvents and hydrogen chloride was passed. In both cases, there 
was no reaction and the ceric oxide remained unchanged. 


1. Reaction in ether—Hydrated ceric oxide (5 g.) and ether (150 c.c.) 
were taken in a flask and hydrogen chloride was passed at 10° C. for three 
hours, with stirring. The reaction was very slow, most of the gas being un- 
absorbed. Finally, the colour of the solution became orange. The flask 
was cooled at 2-3°C. The colour of the solution slowly became deeper 
and at the end of a week, about 5 c.c. of an oily red liquid got collected at the 
bottom. The ether layer had become red in colour. After another week, 
more of the red liquid had formed and then the contents of the flask were 
filtered. The residue on the filter was found to be mostly unreacted ceric 
oxide. The clear filtrate consisting of about 10c.c. of the red oily lower 
layer and about 120 c.c. of the upper layer, was kept at 0° C. for several days. 
Crystallisation did not take place. All attempts to crystallise the substance 
were unsuccessful. 


The two layers were separated and portions of each analysed for cerium 
and active chlorine. Both layers gave the atomic ratio of cerium to active 
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chlorine as 1: 1. It was clear that small amounts of the hexachloroceric 
acid were formed in this solvent. 


A portion of the lower layer was added to some dioxane. The two 
liquids mixed well on shaking. The solution was then kept at 2-3°C., 
when orange red, needle-shaped crystals separated out. The crystals were 
removed, washed with petroleum ether, dried and analysed. The composi- 
tion of the crystals was found to agree with the formula H,CeCl,. 4 C,H,O,. 
Hence it can be assumed that the red liquid that separated contained hexa- 
chloroceric acid, probably with ether of solvation. 


2. Reaction in ethylene glycol mono-methyl ether.—When this solvent was 
used as the reaction medium in the above reaction, more vigorous reaction — 
and quicker changes were noticed. Within about half an hour, nearly 
the whole of the ceric oxide dissolved and the solution became red and 
viscous. The temperature of the flask rose to about 60°C. The flask was 
disconnected from the hydrogen chloride system and the contents filtered. 
There was very little residue left on the filter and a clear, red viscous 
solution was obtained as the filtrate. Analysis showed that the gram 
atom ratio of cerium to active chlorine was 1: 1. 


Attempts were made to isolate a pure compound from this solution. 
Keeping the solution at 0° C. for several days did not effect crystallisation. 
Cooling to — 15°C. also did not make any change. Further cooling in 
solid carbon dioxide bath caused solidification of the entire solution without 
any separation of the acid. The solution was washed with petroleum ether. 
A larger number of washings with petroleum ether gave a pasty product. 
A pure product with a constant composition could not be obtained. 
Removal of solvent by evacuation was also not successful in getting the pure 
acid. 

A portion of the solution was mixed with an equal volume of dioxane 
and kept at 2-3°C. Within a few hours needle-shaped orange coloured 
crystals appeared, which were removed, washed with petroleum ether, dried 
and analysed. The results indicated that the compound had the formula 
H,CeCl, .4C,H;O,. These facts clearly show that hexachloroceric acid 
is formed during the above reaction in presence of ethylene glycol mono- 
methyl ether but the solid acid could not be crystallised. 


SUMMARY 


(1) Benzene, carbon tetrachloride, formamide, cyclohexanol, anisole and 
aceione are not suitable media for the reaction of hydrated ceric oxide and 
hydrogen chloride. 
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(2) The reaction is found to take place to a limited extent in ether and 
nearly to completion in ethylene glycol mono-methyl ether forming hexa- 
chloroceric acid. 
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1. INTRODUCTION 


As early as 1942, Buerger (1942) has pointed out that the Difference- 
Patterson or the Difference-implication diagram could be used to find the 
positions of atoms in structure. Donnay and Buerger (1950) have computed 
the Difference-implication map to get the positions of Mg and Fe atoms in 
the structure of tourmaline. An extension of the Difference-Patterson 
(D-P) was applied by Frueh (1953) to facilitate the solution of order-disorder 
problems in minerals. Again, Buerger (1953) has suggested that an applica- 
tion of the minimum function method to the D-P maps could help in solving 
the crystal structure. 


In a recent paper, Kartha and Ramachandran (1955) have shown that 
the D-P method can be systematically and successfully applied to analyse 
crystal structures. The present paper is mainly concerned with an attempt 
to apply the method to a more complicated example, namely, the isomorphous 
crystals di-p-tolyl telluride and di-p-tolyl selenide, whose structures have 
been worked out by Blackmore and Abrahams (1955 a, 5). In their report, 
complete data are given of all the measured structure factors and these could 
be utilised to calculate the D-P diagram. The two compounds -crystallise 
in the orthorhombic system, space group P2,2,2,, with four molecules per 
unit cell. The unit cell dimensions are: 


Selenide: a= 25-12A; b=7:99A; c=5-88A. 
Telluride: a= 25-33A; b=8-05A; c=6-01A. 
2. THe D-P DIAGRAM OF THE b PROJECTION 


The electron density projections given by them showed that the b-axis 
projection was free of overlap and so it was first tried. Using the data of 
F (h0/), |F (A0/) |? were calculated for the two crystals and the Patterson 
function was computed in each case. The Patterson function of the selenide 
was subtracted from that of the telluride and the difference was plotted. 
The resulting D-P diagram is shown in Fig. 1. 
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Fic. 1. The D-P projection of di-p-tolyl telluride and selenide along the b-axis. Contours 
for the heavy atom peaks are at intervals five times as much as the intervals of contours for other 


peaks. 


The D-P diagram should contain only peaks at the termini of vectors 
joining the replaceable atoms (Te or Se) with themselves or the other atoms 
in the structure. The strong peaks, corresponding to the vectors joining 
the replaceable atoms only, occur in the diagram at the positions: (1) 0,0; 
(2) 0:25, 0-30; (3) 0-25, 0-50; (4) 0-25, 0-70; (5) 0-50, 0-20; (6) 0-50, 
0-80; (7) 0-75, 0:30; (8) 0:75, 0:50; (9) 0-75, 0-70. Of these, peaks 
2,4, 7 and 9 are of nearly the same strength while 3, 5, 6 and 8 are of double 
the strength, while peak 1 at the origin is four times as strong. The positions 
of the replaceable atoms could be obtained from the D-P diagram by trans- 
lating it so as to make one of the peaks (e.g., 2) coincide with the origin and 
marking out the positions of the peaks which superpose (the theory of this 
is simple). Since the structure has a centre of symmetry in the projection, 
this straightaway gives the positions of the replaceable atoms, namely 1, 3, 5 
and 9. These are designated by a,, dz, a; and a, (Fig. 1). 


In order to get the positions of the other atoms, the D-P diagram was 
then given a translation equal to (a, 4,), (a3 a,) and (a, so as to bring 
each of the other peaks in coincidence with the peak a, and the positions of 
the peaks which superpose were marked in each case. Owing to diffraction 
effects some of the fainter peaks were suppressed or displaced if they occurred 
near a strong peak in the D-P diagram. Consequently an average was taken 
of all the measurements. The molecule thus determined is shown associated 
with each of the replaceable atoms in Fig. 1. 


In finding the superposition it must be remembered that in the D-P 
diagram a certain peak will be doubled or quadrupled, according to the 
teplaceable atom peak with which it is associated. A certain amount of 
judgment is therefore needed in finding the superposed positions. 
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3. COMPARISON WITH EXACT STRUCTURE 


The fractional co-ordinates of the atoms in the molecule with respect 
to the heavy atom are listed in Table I, together with the exact values for 


TABLE | 
D-P method Electron density method 
Atoms Telluride Selenide 

x x’ x 2” 
Cc, —0-060 0-217 —0:0614 0-2161 —0:0604 0-1918 
Cc, —0-056 0-408 —0:0545  0-4133 —0:0551 0-3972 
—0-100 0-600 —0:0961 0-5602 —0:0994 0-5367 
Cc, —0-148 0-500 —0-1458 0-5144 —0-1484 00-4700 
C; —0-150 0-292 —0-1528 0-3160 —0-1539 0-2649 
Cc —0-100 0-175 —0-1108 0-1691 —0-1097 00-1262 
C, —0-190 0-733 —0-1913 0-6733 —0-1955 0-6187 
C,’ +0-060 0-217 +0-0637 0-2128 +0-0602 0-2020 
C,’ +0-056 0-408 +0-0637 0-3960 +0-0558 
C,’ +0-100 0-600 +0-1074 0-5428 +0-1002 0-5540 
Cc,’ +0-148 0-500 +0-1503 0-5016 +0-1481 
+0-150 0-292 +0-1503 0-3201 +0-1516  0-2833 
C,’ +0-100 0-175 +0-1074 0-1732 +0-1077 0-1391 
C,’ +0-190 0-733 +0-1975 0-6583 +0-1957 0-6459 


the telluride and selenide structures as determined by Blackmore and 
Abrahams (1955 a, 5). It will be noticed that the agreement is fairly good 
except for the methyl carbon atoms C; and C,’ for which the z co-ordinates 
are appreciably different. This is due to the diffraction effect of the strong 
peaks. In all the eight positions, this particular atom is displaced roughly 
in the same way so that the displacement is not evened out by averaging. 


Another point worthy of note is that the molecule as determined from 
the D-P diagram is symmetric about a line parallel to the c-axis. This sym- 
metry is really absent, but in the superposition method of finding the peaks 
the differences noticed were not considered significant and. an average was 
struck between the two halves of the molecule. 


Further it must be mentioned that the Se and Te atoms are not exactly 
in the same positions in the two crystals. In fact a typical Te-Te vector has 
the components 0-26, 0-28 (as determined from the corresponding simple 
Patterson function) while the corresponding Se-Se vector is 0-26, 0-22 and 
in the D-P diagram the corresponding peak occurs at 0-25, 0-30 which is 
to be expected. So it is interesting to note that the D-P diagram has quite 
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sharp peaks corresponding to the replaceable atom interactions, in spite of 
the relative displacement of the corresponding atoms in the two isomorphous 
crystals and that the weaker peaks corresponding to the other interactions 
are also well defined. This shows that the D-P method can be confidently 
applied in general to isomorphous crystals even though the replaceable atoms 
may slightly differ in their positions. 

4. APPLICATION OF VECTOR SET METHOD TO THE PATTERSON DIAGRAM 


During the preparation of the Difference-Patterson diagram, it was 
observed that there was a close correspondence between the peaks in it and 
the peaks found in the ordinary Patterson diagram of the selenide and the 
telluride. This is to be expected because both Se and Te atoms are so much 
larger scatterers than the carbon atoms that the C-C interaction will be of 
small intensity compared to the Te- or Se-C and Te-Te or Se-Se interactions. 
Consequently the methods described in the previous section’ was applied 
directly to the Patterson projection of di-p-tolyl selenide along the b-axis. 
But for a background due to the C-C peaks which is seen from the fact that 
the minimum value of the Patterson function is nearly 2,000, as compared 
with value close to zero of the D-P map, this diagram is not inferior to the 
D-P diagram. However, there is the advantage that the shift and a small 
broadening of the peaks, which is present in the D-P diagram due to the 
difference in the positions of the atoms in the two crystals, is absent here. 


Table II gives the co-ordinates of the atoms in one molecule obtained 
from it. The deviation 4’ from the exact value are shown in the table and 
may be compared with the deviations 4 found with the structure deduced 
from the D-P. It will be observed that the deviations are small and more 
or less the same but for the deviations in the position of the methyl carbon 
atoms, which are fairly large in both cases. 


In order to test the usefulness of the method, the R value was calculated 
with the co-ordinates in Table II and was found to be 0-36 taking all the 
reflexions observed by Blackmore and Abrahams (1955 5) in this zone. It 
may be mentioned that taking the co-ordinates from the telluride structure 
unchanged the above authors found an R value of 0-45. Thus the simple 
Patterson could have been used directly to get an even better first approxi- 
mation by the application of vector set method. In fact if it is assumed 
that the methyl carbon atom is in line with the two benzene ring carbon 
atoms, which are in the para positions relative to each other, then the mean 
deviation would be even less. 


The same method was applied to the Patterson projection of di-p-tolyl 
selenide along the c-axis. The projection is shown in Fig. 2. As before 
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Fic. 2. The Patterson projection of di-p-tolvl selenide along the c-axis. Contour intervals 
areas in Fig. 1. 


the strong peaks are designated by numbers from | to 7, and the positions 
of the replaceable atoms by 4), a2, ad, and a,. The y co-ordinate thus obtained 
is very close to 0-25. Because the y co-ordinate is close to 4 every replaceable 
atom peak is doubled and further the rest of the molecule associated with 
the two peaks are related by a reflection on the b-plane. Hence it is difficult 
to uniquely interpret the diagram. The two possible positions of the molecule 
associated with one of the strong peaks are shown in figure. However if 
the presence of a benzene ring is previously known, a rough idea of the dis- 
position of the molecule can be got. This projection is therefore more diffi- 
cult to interpret than the previous one. It is thus seen that the interpretation 
of the heavy atom diagram or the D-P diagram may become difficult if the 
heavy atom or replaceable atom has special values for the parameters. 


Thus, in general, it is clear that Buerger’s idea of solving the Patterson 
map directly for the crystal structures can be profitably applied if there is a 
heavy atom and if there is no appreciable overlap in the projection. In 
the present case the structure had a centre of symmetry in the projection and 


no doubling of structure was introduced. If such a symmetry is not present, 


then the structure can be uniquely solved if (a) the number of replaceable 
atom peaks is greater than two, and (b) the group of replaceable atoms do 
not have a centre of symmetry (Kartha and Ramachandran, 1955). 


The D-P technique is being systematically applied for solving unknown 
crystal structures in this laboratory. 


In conclusion, the author wishes to express his sincere gratitude to 
Professor G. N. Ramachandran, for the keen interest he took in this study. 
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SUMMARY 


The Difference-Patterson (D-P) technique is applied to a complicated 
example, namely, the isomorphous crystals, di-p-tolyl telluride and di-p-tolyl 
selenide, whose structures have been recently reported. It is shown that the 
D-P method can be confidently applied in general to isomorphous crystals, 
even though the replaceable atoms may slightly differ in their positions. 
Also, Buerger’s vector set method is applied to the Patterson projections of 
the di-p-tolyl selenide along b- and c-axes. It is shown that Buerger’s method 
can be profitably used when there is no appreciable overlap in the projection 
and if the heavy atoms or the replaceable atoms do not have special values 
for their parameters. 
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An adduct, formed from f-naphthol and maleic anhydride, by heating an 
equimolecular mixture at 220-40° C., was first isolated by Takeda and 
Kitabonoki.! By Clemmenson reduction of the adduct to a desoxy-com- 
pound, identical with the product of hydrogenation of the adduct of naphtha- 
lene and maleic anhydride, they showed that the structure of this adduct 


was represented by I. 


Later Takeda, Nagakura and Kitabonoki? isolated another adduct 
from the reaction product obtained by heating f-naphthol and maleic an- 
hydride at 220-40°C. The steric configurations of these two adducts were 
then proposed on the basis of their dipole moments; the one with the smaller 
moment (2-86 D) being assigned the exo-configuration (II) and the one with 
the larger moment (6-54 D) the endo-configuration (III). They argued that 
in the exo-form the dipole moment of the carbonyl group (2-5 D) and the 
acid anhydride (4-32-4-05 D) will be in opposite directions while in the endo- 
form these two would be in the same direction. 


endo 
The large difference in the two dipole moments and the good agreement 
between the two observed and calculated values makes the correctness of these 
assignments likely, but the misleading deductions possible from dipole moment 
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measurements in such ring systems * * makes a rigorous chemical proof 
desirable. This has now been provided by the following transformations. 
The two anhydrides were converted into the acids by boiling with water and 
then into the respective esters by diazomethane. These isomeric keto esters 
were then reduced by sodium borohydride in methanol. The endo-form gave 
a lactone ester C,;H,,0, (V) by elimination of methanol, while the exo-form 
gave the normal reduction product C,,H,,O; (VI). The endo acid obtained 
by boiling the endo adduct with water was also reduced with potassium 
borohydride in aqueous medium and the resulting lactone acid (IV) converted 
into its ester by diazomethane. It was found that this ester was identical 
with the lactone ester obtained by the reduction of the dimethyl ester of the 
cndo acid. The scheme is given below. 
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The peak inthe infra-red spectrum of the lactone ester (Fig. 1) at 
1776 cm.— shows the presence of the y-lactone ring. The formation of this 
lactone ester proves conclusively the structure of the endo adduct. 


EXPERIMENTAL 


Preparation of adducts 

Purified 8-naphthol (9 gm.) and pure maleic anhydride (6 gm.) were 
heated at 220—-40° C. in a conical flask placed in a Wood’s alloy metal-bath 
for half an hour. After cooling, the products were dissolved in the minimum 
quantity of hot ethyl acetate and allowed to crystallise slowly. Crystals 
obtained melted at 152-58° C. (4-5 gm.). On recrystalling thrice from ethyl 
acetate, pure crystals of the exo-form (1-5 gm.) melting sharply at 193-94° C. 
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were obtained. Careful recrystallisations of the residues from the above 
yielded the endo-form (0-4 gm.) melting at 192~93° C. (mixed melting point 
with the exo-form 152-55° C.). 
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Incidentally it was found that the mixed anhydrides could be converted 
into the dicarboxylic acids by heating with water and separated into the exo 
acid and endo acid easily. The former, unlike the latter, was only slightly 
soluble in ether. 


Conversion of the adducts into the dimethyl ester 


The dicarboxylic acids (0-5 gm. each) were dissolved in methyl alcohol 
(Sc.c.) and an ethereal solution of diazomethane was added until a yellow 
colour persisted. The yellow colour was then destroyed with a few drops of 
glacial acetic acid. On concentrating the solutions and cooling, crystals of 
the dimethyl ester separated (0-40 gm.). On recrystallisation from methanol, 
the exo-ester melted at 96-97° C. and the endo-ester at 112-13°C. The ultra- 
violet spectrum of the anhydrides (Figs. 1 and 2)—A,,.,, of exo=297 mp and 
e = 305-8 and A,,.,, of endo=299 mp and « = 338 as well as the melting 
points of the esters show indentity with the values given by Takeda et al. 
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The anomalous, long wave-length band in the U.V. spectrum due to the inter- 


action between the carbonyl group and the benzene ring will be discussed 
elsewhere. 


Reduction of the esters 


(A) Exo-ester.—The ester (0-70 gm.) was dissolved in methanol (10 ml.) 
and sodium borohydride (0-5gm.) was added in small quantities. A 
vigorous reaction took place. After two hours, water was added and then 
extracted with chloroform. On removal of the chloroform, an oil (-50 gm.) 
was obtained. The oil was dissolved in 95 per cent. ethanol and cooled. 
Colourless crystals melting at 118-19° C. were obtained. Recrystallisa- 
tion from alcohol gave needles melting at 128-29°C. C,.H,sO; requires 
C = 66-20 per cent.; H=6-21 per cent. Found: C = 66-32%, H = 6-24%. 


(B) Endo-ester.—The ester (0-75 gm.) was dissolved in methanol (10 ml.) 
and worked up exactly as in the above case. Pure crystals (0-25 gm.) melting 
at 128-29°C. were obtained. Mixed melting point with the reduction 
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product of the exo-ester was 95-97°C. (C,;H,,O, requires C = 69-75%, 
H =5-46%. Found: C = 69-76%, H = 5-37%. 


Reduction of the endo-dicarboxylic acid 


The acid (1-0 gm.) was dissolved in aqueous methanol and potassium 
borohydride (2-0 gm.) was added in small quantities. A vigorous reaction 
followed. After two hours, the resulting solution was just acidified with 
glacial acetic acid. The solid, which separated, was filtered out and the 
filtrate heated with a few drops of concentrated hydrochloric acid. On 
allowing the solution to cool, needle-shaped crystals (0-3 gm.) separated. 
Recrystallised from aqueous methanol the lactone acid melted at 198-200° C. 
C,,H,,O0, requires C = 68-85%, H = 4-:92%. Found: C = 68-88%, H = 
4-97%. 


These crystals (0-2 gm.) were dissolved in methanol and esterified with 
diazomethane in ether in the usual way. On-removing most of the solvents 
and leaving overnight, crystals melting at 114-16°C. (0-10 gm.) were 
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obtained. After recrystallisation from methanol, these melted at 128° C,, 
and gave no depression in melting point with the lactone ester obtained from 
the reduction of the endo-dimethyl ester by sodium borohydride. 


The author wishes to record his deep indebtedness to Dr. R. C. Cookson, 
Senior Lecturer in Chemistry, Birkbeck College, London, for valuable 
suggestions during the course of this work. His thanks are also due to 


the University of Travancore, for the study leave granted to him which 
enabled him to do this work. 
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(Communicated by Prof. B. R. Seth) 
1. INTRODUCTION 


Tue flow of a viscous imcompressible fluid in an elastic tube has been studied 
by Rashevsky (1945) and Morgan (1952). Morgan has examined the condi- 
tions under which an approximate relation between the radius of the tube 
and the distance along the axis is consistent with the assumptions made in 
the solution. The solution obtained is limited to Newtonian liquids and 
do not give any information about the effect of cross-viscosity which is 
present in a number of highly viscous fluids of technological and biomathe- 
matical importance. 


In this paper the flow of certain non-Newtonian liquids through an 
elastic tube is studied by considering the general Reiner-Rivlin stress-strain- 
velocity isotropic relations instead of stress-strain relations of classical 
Hydrodynamics. The most general form of such relations is obtained by 
introducing second degree terms of the type d,‘d;* into the classical relations. 
The conditions are examined under which an approximate relation between 
the radius of the tube and the distance along the axis is consistent with the 
assumptions madé in the solution. These conditions are reduced to the 
relation between three dimensionless parameters. Second approximations 
are considered for different ranges of values of these parameters. It is found 
that the effect of cross-viscosity is to increase the radius of the tube with 
the distance along the axis. 


2. EQUATIONS OF MOTION 


The equations of motion proposed by Reiner (1952) for non-Newtonian 
liquids are 


dv; 
p + 04,524) (2-1) 


0, (2-3) 
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where f* is the body force, 2dj'= v;,;-+-j,; is the strain velocity tensor, »; js 
the velocity vector and p is the density (assumed constant), F, and F, are 
the two constants, known as the coefficients of viscosity and cross-viscosity 
respectively. 

For the problem on hand it is convenient to use the cylindrical co-ordinate 
system r, 9, z in which r is the distance from the axis of the tube and z is 
the distance along the axis. The velocity components at any point are 
given by v;, vg and vz. On account of axial symmetry, the velocity com- 
ponent v, and its derivatives vanish. 


The equation of continuity is: 


2 (rvy) 
r 


+ = 0. (2.4) 


The equation (2-4) can also be written in an integral form, viz., that the 
total mass flow across any cross-section of radius R (z) is constant and 
equal to say, Q: 


R(zd 


rvzdr = Q. (2.5) 


The equations of motion for axially symmetric, steady laminar flow of a 
non-Newtonian liquid become 
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where p is the pressure (for convenience p is taken as the difference between 
the fluid pressure and the constant pressure outside the tube), » = p/p, 
», =}4/p are the Kinematic coefficients of viscosity and cross-viscosity res- 
pectively. On using (2.2) the normal compressive force per unit area per- 
pendicular to a radius is given by 


Ere 2) =p — 2u — [ + + 2-8) 


EQUATION OF ELASTICITY 


If R varies very slowly with z and if we neglect effects at the ends of the 
tube, then the forces in the tube will have no component parallel to the axis. 
The following equation can be derived, provided 5 the wall thickness, is 
very small compared to R: 


(2-9) 


where E is Young’s modulus, Ry, is the unstretched radius of the tube (subse- 
quently assumed to be constant) and E,, represents the difference between 
the normal compressive stress on the interior of the tube and the external 
pressure. 


In this derivation we have neglected the effect on the radius R (z) of the 
fluid friction on the interior of the tube, and that, due to the slow variation 
of R(z), we take the component of E,, which is normal to the wall to be 
equal to E>, itself. 


3. SOLUTION OF THE EQUATIONS OF MOTION 


First approximation.—Since R (z) varies very slowly with z, we can to 
a first approximation neglect all terms containing v, and its derivatives. 
We will investigate under what conditions this assumption may be expected 
to yield reasonable results. Equation (2-6) becomes 
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Now neglecting (=) we get 


i.e., p=p(z). 
Hence from (2-8) it follows that 

Err (R, 2) = p(2). 
Then (2-9) gives 


dp __ EsdR 
dz R®dz’ 


and equation (2-7) can now be written 
1 dw, 
= Tr sf); 


provided we neglect terms 


dvz 


dvz, Le dvz 
droz’ roo oz? 


The integration of (3-2) and the application of the boundary condition 
vz = 0 at r=R yields 


(3.3) 
For any fixed z this equation gives a parabolic radial distribution for vz as 
in the case of classical viscous fluid. 
Substituting (3.1) into (3.3) we obtain 
r2 
If we now apply equation (2-5), we obtain 


dR _ — &Q 
dz ~ EdapR® 


_ 


_ 
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as was obtained by Morgan (1952) and Rashevsky (1945) in the case of 
classical viscous fluid, R, being the radius at z = 0. 


Validity of the approximation—We have neglected certain terms in 
equation (2-7). We will now examine under what circumstances the above 
solution is consistent with the assumptions. 


Now from equations — (3-2) and (3-5) we have 


row mp*R** 


=0O aR) ’ 
3? _ 


1 dvz _ pp, Q3 ) 


Hence We may neglect ? and “2 provided 


1. (3-7) 


If we may write Q = pVR* where V is the average velocity over a cross+ 
section, we have 


(3-8) 
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< 1. 


We may neglect the terms 


r ow oz’ 


<1, 


v2 
ER <1. (3.12) 


We may neglect the term », we. = provided, 


<1 (3-13) 


<!- (3.14) 


Now if we write pvr = R,, a Reynolds number appropriate to our 


problem (V and R may be taken to be the values at z = 0) and 

p4V? 

ES ESR 

the conditions (3-8), (3-10), (3-12) and (3-14) reduce respectively to 
< 1, R,? < 1, Rs < 1, 1. 

But these reduce to only three conditions 


RR, <1, (3-15) 
l, (3-16) 
R; <1. (3-17) 


R, and are non-dimensional parameters. R, = ratio of a viscous to 
an elastic force. &; = ratio of a cross-viscous to an elastic force. We 


regard R,, R, and R, as the three dimensionless parameters of the system. 
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For any fixed 2%, the flow will approach Poiseuille flow as ®, and 
0. 
Now considering equation (2-8) we can write 


| | < | az |’ 


| 
| S| 32 


and 


ie., provided <1, R2R; <1 and <1, which have already been 


obtained as conditions. Hence we are justified in writing 
in the case of classical viscous fluid. 
Using equation (2-6), we have, 


20Q, + higher order terms. 


By integrating from r to R and using (2-9) we get 
Fr 
P(r,2) = ™ R): G-18) 
By differentiating this equation (3-18) with respect to z we get 


Thus the expression for 2 will be the same as was obtained by Morgan 


(1952) in the case of classical viscous liquid. The last two terms in equation 
(3-19) are negligible in comparison with the first term provided R,*< 1, 
which is the condition already obtained. 


Second approximation—We will consider three cases for ranges of 
values of the three dimensionless parameters of the system. 


Cases : 1, RRy<K<1, R <1. 
ad R= O(1), 
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In these cases we will have to take account of some of the terms pre- 
viously neglected, by evaluating them with the help of the first approxima- 
tion. Calling the terms retained in the first approximation, terms of order 
one, the next higher order terms in these cases will be those of order &,&,, 
ie., the quadratic terms in (2-7). These terms have the same approximate 
value as in the case of classical viscous liquid. That is, 


d 2 dR 
Hence (2-7) in this second approximation becomes 
1 2 d2u 12 


Integrating this and applying uh boundary condition vz= 0 at r= R, we 
obtain as Morgan (1952) has done in the case of classical viscous liquid, 
Now applying (2-5) we arrive at the same equation for a as was obtained 
by Morgan (1952) in the case of eet viscous liquid, viz., 
Edzp 


Since = is negative, the effect of the second term in parenthesis is to 


diminish the radius R more rapidly. It is to be noted that, since the condi- 
tions (3-7), (3-9) and (3-11) involve the variable R the equations hold only 
over a range of z in which R does not become too small. 
Case: R, <1, R= O(R,) and R = O(R)). 

In this case we must take account of terms of order R,R,, R,? and Rs, 
i.e., all the terms previously neglected in (2-7) and (2-8). We will follow, 
as before, the method of successive approximations, i.e., by replacing all 
unknown terms of second order by their approximate values as obtained 
from the first approximation. 

The equation = 7) becomes 

8Q? dR , 1 dy 

32pr,Q? dR 
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Now equations (2-6), (2-8) and (2-9) yield the following expression for p: 


4 a (a 


3QESu, dR 
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dz?\R# (3-26) 


If we introduce this expression (3-26) for P into (3-24) and perform 
the integration, applying the boundary condition, viz., when r = R, vz = 0, 
then we get by expressing = and 7" r in terms of = A 
+ de (— Re + 


mit “dz (RO 


(3-27) 
Now the application of equation (2-5) will yield the equation 


Eézp Q , 52 1 4u,Q? | 


(3-28) 
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The last two terms in parenthesis in equation (3-28) are the cross- 
viscous terms. While the effect of inertial forces and viscosity is to dimi- 
nish the radius of the tube, the effect of cross-viscosity is to increase the radius 
of the tube with the distance along the axis. Thus it is seen that the effect 
of cross-viscosity compensates to some extent the effect of viscosity and 
inertial forces in flow of a liquid through an elastic tube. Prof. G. W. 
Morgan has kindly pointed out to me that the influence of Poisson’s ratio 
on the shape of the tube might be as important as the effect of cross-viscosity, 
Since the purpose of our problem is to investigate only the effect of cross- 
viscosity on the shape of the elastic tube this effect has not been taken into 
consideration. 

SUMMARY 

This paper is a study of the effect of cross-viscosity on the radius of an 
elastic tube when certain highly viscous non-Newtonian liquids flow through 
it under pressure gradient. The effect is found to be of opposite nature to 
that of inertial and viscous forces and causes the radius of the tube to in- 
crease with the distance along the axis. 
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(Communicated by Sir C. V. Raman) 
1. INTRODUCTION 


THE treatment customarily adopted’ for analysing the propagation of light 
in optically active crystals is somewhat lengthy and cumbersome, and is even 
so not free from approximation: its mathematical complexity being perhaps 
partly the reason why an alternative quasi-theoretical analysis not involving 
the electromagnetic theory is sometimes used.*»* A concise yet logical 
treatment of the problem is, however, easily possible without the introduc- 
tion of any fresh physical hypothesis—as will be shown in this paper. 


We have merely to remember that for an electromagnetic wave inside 
a crystal, it is the dielectric displacement vector D (the ‘ vibration’), and not 
the electric intensity E which necessarily lies on the wave-front. Hence 
we seek to express E in terms of D, rather than vice-versa. Even for the case 
of optically inactive crystals, such a procedure—though not commonly 
adopted—is known to simplify the theoretical treatment, as may be seen 
especially in the case of absorbing crystals‘; and the geometric representation 
of the optical properties of a crystal by an index ellipsoid follows directly 
in the structure of such a presentation instead of having to be proved by 
indirect means later—as is usually done.® 


2. FORMULATION OF THE PROBLEM 
Any wave travelling in the crystal along an arbitrary direction ‘ z’ (con- 
veniently taken as being normal to the plane of the paper) may be specified 
by giving the x and y components of the D-vibration: Dz exp iw (t — z/v) 
and Dy exp iw(t — z/v). Here v is the velocity of the wave and the ratio 
(D,/Dz) completely defines the polarisation of the wave. Owing to the re- 
quirement that the vectors D and E must satisfy both Maxwell’s relations 
as well as the characteristic polarisable property of the medium, only two 
specific waves can be propagated in the z-direction—the problem being to 
determine their velocities, v, and v,, and their states of polarisation 
(D,/D;), and (D,/Dz)s. 
247 
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Maxwell’s equations for a plane-wave field reduce to the following 
single equation when H is eliminated*—as is shown in standard texts’ 


= c? [E — s(Es)] (1) 
where s is a unit-vector along the wave-normal. For inactive crystals the 
polarisable characteristics of the medium lead to the usual relation D = (¢) E, 
where (¢) represents the symmetric dielectric tensor—which operates on E 
to give D. But for optically active media it is assumed that 

D=(«e) E+iGxE (2) 
where the additional antisymmetric term represents a feeble additional 
polarisation which oscillates out of phase with the electric intensity (as indi- 
cated by the factor i), and yet does not introduce any absorption because it 
is also orthogonal to the electric intensity (as is indicated by the occurrence 
of the vector product). The gyration-vector G corresponding to the direction 
of propagation s determines the optical rotatory power in that direction. 


3. THE INVERSE GYRATION- VECTOR 


It is known that by taking the z-axis along the wave-normal (instead of 
choosing co-ordinate axes along the electrical axes, as in usual treatments), 
we not only confine ourselves to two components of D instead of three; in 
addition, the components of the vector equation (1) take the following intelli- 
gible form: 


= c*Ez; v*Dy = c*Ey; Dz =0 (1 a) 
where, it may be noted, separate and explicit expressions are obtained for 
Ex and Ey. 

Equation (2) which expresses D as a vector function of E, is now trans- 
formed by expressing E in terms of D using the inverse vector function. 
Since this inverse function will obviously contain not only a real symmetric 
part but also an imaginary antisymmetric part (giving that part of the electric 
intensity which oscillates out of phase with D, and is orthogonal to it), we 
may write 

=(a)D—ilxD (2 a) 


The first term alone would be present for an optically inactive crystal, the 
symmetric tensor (a) then defining the index ellipsoidt—being proportional 
to the inverse of the dielectric tensor (see Appendix). The comparatively 


* We are not here considering other theories of optical activity?» in which B + H. 


+ If ayyx? + aooy? + ....2a3, zx = 1, be the equation to the index ellipsoid, then the a,; are 
also the components of the tensor (a) with respect to the co-ordinate axes chosen. 
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small second term (which introduces the optical activity) is determined by 
the new vector I"—which we shall refer to as the inverse gyration-vector 
corresponding to the direction of propagation s. 


4. DERIVATION OF THE EQUATIONS OF WAVE-PROPAGATION 


Writing down the x and y components of equation (2a) and omitting 
the terms containing the factor Dz (since Dz = 0): 


= ayDz + + 
c*Ey = y + 42D z = 


Here a,x” + doay? + 2a;,xy = 1, is the equation to the elliptical section 
of the index ellipsoid taken perpendicular to the wave-normal. We now 
choose the x and y axes parallel to the principal radii of this elliptical section, 
which represent the vibration-directions in the absence of optical activity. 
With this understanding we may, in Equation (3), set dy, = 0, a,, =v’ and 
dg, = v''*; here v’ and v” represent the velocities of the waves in the absence 
of optical activity, because of a well-known property of the index ellipsoid 
which may also be obtained by setting I’, =0 at the end of the present dis- 
cussion. Substituting for c*Ez and c?Ey from equation (1 a) we thus obtain 
as our fundamental equations: 


(3) 


(4) 


These fundamental equations are essentially the same as those obtained in 
Pockel’s Lehrbuch* starting from an entirely different theory of optical 
activity. There are two pairs of roots v,, (Dy/Dz), and v,, (Dy/Dz)., which 
simultaneously satisfy these equations—giving thereby the required velo- 
cities and states of polarisation. 


5. THE VELOCITIES AND STATES OF POLARISATION OF THE WAVES 


Multiplying the two equations of (4) to eliminate (D,/Dz) we obtain the 
following quadratic in v*, whose roots v,* and v,? give the velocities of the 
waves : 


(v? — v') — =? (5) 
Subtracting the second equation of (4) from the first to eliminate v?, 
we get: 
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This is a quadratic in (D,/D), and its roots uniquely determine the states of 
polarisation because of the relation: 


Dy _ |Dy! 


where the ratio | Dy | : | Dz | is obviously the ratio tan @ of the (real) ampli- 
tudes of the y and x components of the D-vibration; and R is the relative 
phase difference between these components. Since the right-hand side of 
(6) is entirely imaginary it follows that both roots of the equation are imagi- 
nary (R= +7/2). Further since the roots of (6) are obviously reciprocals 
of one another, it follows that if one root is of the form (B,/A,) e**’®, the other 
will be (A,/B,) e-**’*. This signifies that both vibrations are lliptically 
polarised with their respective major axes along the two perpendicular princi- 
pal planes, the vibrations being of the same ellipticity but traced in opposite 
senses. 


The value of tan 24, which determines the common ellipticity of the two 
vibrations is obtained directly from equation (6) by substituting in it the 
value of (Dy/D.) which is equal to (+ i tan 9,) for the right elliptic vibration 
and (— i tan for the left elliptic vibration: 


| tan 20 | = — | 


(7) 


Similarly the difference in the squares of the velocities of the two waves 
may be obtained from (5)—-since it is the square-root of the discriminant: 


(v1? — 07) = — + (202)? (8) 


In the absence of linear birefringence (i.e., if we imagine v’ = v”), the 
waves will be circularly polarised according to (7), the rotatory power p being 
obtained from (8) as: 


p = (w/vm*) - Tz (9) 

where Um is a mean velocity (see Pockels*). 
We may mention that because of an approximation made in the usual 
treatment, its resultant equations differ slightly from (5) and (6); also, equa- 


tions (7) and (8) are the accurate forms of the approximate relations obtained 
by the method of superposition.*»* 


: 
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6. THE VARIATION OF ROTATORY POWER WITH DIRECTION 


The gyration vector G, and hence in the present treatment the inverse 
gyration-vector J", are both vectors not having their axes fixed with respect 
to the crystallographic axes,{ but are linear vector functions of the direction 
of propagation. Thus G = (g)s and I’=(y)s, where the gyration tensor 
(g) and hence the modified gyration tensor (y) are both nine-component 
tensors. The rotatory power is however determined, according to (9), only 
by the resolved component I’, of the inverse gyration-vector along the direc- 
tion of propagation. (We are here denoting this by I’; instead of Iz since 
we are choosing an arbitrary co-ordinate system x’y’z’.) Since an analogous 
result holds in the usual treatment it follows as in that case® that this scalar 
parameter of rotation I’; is a quadratic function of the direction cosines 
I'm'n' of propagation: 


Ps = yyyl'?+ + 2yegm'n' + 
where 
ij + Vii)» (i, i, 2, 3). 


The discussion of optical activity in relation to crystal symmetry follows 
as in the usual treatment. 


Finally, the author wishes to express his gratitude to Prof. Sir C. V. 
Raman, F.R.S., N.L., for his kind interest in this investigation. 


7. SUMMARY 


In the usual treatment of the optical activity of crystals, the displacement 
vector of the light wave is expressed as a function of the electric vector—this 
being done with the aid of the dielectric tensor and the gyration-vector. The 
treatment is made more concise and tractable by using the inverse vector 
function and expressing the electric vector in terms of the displacement 
vector, since it is the latter which necessarily lies on the wave-front; this 
can be done by using the inverse of the dielectric tensor (which defines the 
index ellipsoid), and an inverse gyration-vector (which determines the rota- 
tory power for the particular direction of propagation considered). The 
inverse gyration-vector (like the gyration-vector) is itself a function of the 
direction of propagation, being related to it by a modified gyration tensor— 
which consequently determines the rotatory power for all directions. 


t The new treatment given in Sommerfeld’s Optics® is deficient in this respect. 
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APPENDIX 


Since the index ellipsoid defined by the tensor (a), and the inverse gyra- 
tion-vector I” may be directly taken as the phenomenological quantities 


required to describe the propagation along any particular direction s, the 
relation between I and the gyration-vector G would not usually be required. 
The relation may however be obtained by writing (2) as D=D’ + iP’, where 


E; = GxE (11) 

Substituting for D in (2 a) 
+ i@P’— xD’ 
Equating real parts and neglecting "x P” in comparison with c*E: 
(a) = (<)> (12) 
Equating imaginary parts: (a)P” =I°xD’. Choosing axes of co-ordinates 
X, Y, Z, along the principal electrical axes of the crystal, the X-component 
of this relation is 
ayP,” = — 
or, by virtue of (11), if «x, €y, €z, be the principal dielectric constants, 
ax (GyEz — G,Ey) = («zEz) — Iz (eyEy) 


We may separately equate the coefficients of Ey and Ez occurring on both 
sides since this relation has to hold not just for one value of the ratio Ey/E;, 
there being two polarised waves that can be propagated along the same 
direction. Hence, using (12), 


2 2 2 


€z€x €x €y €y €z 


The relation between the components of the modified gyration tensor and. 
those of the gyration tensor may now be easily written down: 


c c? 
Vai = Yu = Sais Yai = Sai (i= 1, 2, 3). 


€x ¢y 
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